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ABSTRACT 



A new classification system for carbon-rich stars is presented based on an 



analysis o: 



available (IKeenan 



51 AGB carbon stars through t he most relevant classifying indices 



1993 



Morgan et al. 



20041 ) . The extension incorporated, that 



also represents the major advantage of this new system, is the combination of 
the usual optical indices that describe the photospheres of the objects, with new 
infrared ones, which allow an interpretation of the circumstellar environment of 
the carbon-rich stars. This new system is presented with the usual spectral sub- 
classes and C2-, j-, MS- and temperature indices, and also with the new SiC- 
(SiC/CA. abundance estimation) and r- (opacity) indices. The values for the 
infrared indices were carried out through a Monte Carlo simulation of the radia- 
tive transfer in the circumstellar envelopes of the stars. The full set of indices, 
when applied to our sample, resulted in a more efficient system of classification, 
since an examination in a wide spectral range allows us to obtain a complete 
scenario for carbon stars. 



Subject headings: catalogs — circumstellar matter — infrared: stars — stars: 
atmospheres — stars: carbon 
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Introduction 



The asymptotic giant branch stars with a ratio C/O > 1 have their optical spectra 
ruled by bands of carbon compounds, which obscure many atomic features. The green-red 
optical spectrum is dominated by Swan bands 12 C 12 C, Red System bands 12 C 14 N and 
sometimes present isotopic bands, e.g. 13 C 12 C and 13 C 14 N. As a result, a classification 
carried out through the optical atomic data on carbon stars is troublesome. During their 
ascent on the AGB phase, the mass loss from the star creates a circumstellar envelope of 
gas and dust. The compounds of this shell have their maximum emission on the infrared 
spectral region. The main infrared feature for carbon stars is the 11.3/zm emission due 
to the presence of SiC grains, which also represent an evidence of a C-rich dust envelope. 
Spitzer IRS spectra also confirmed two absorptions at 7 and 13.7 fim, whose origin is still 



not well understood (ISpeck et al. 



2006|). 



The first classification system for carbon-rich stars was the Henry Draper Catalogue 



(jCannon &: Pickering 



19181 ). in which the stars were presented in two spectral classes: type 
R and N, that were also divided into temperature subclasses. These subclasses were based 
on lines in the blue spectral region, that, although being a good source of information to 
determine the temperature classes of G and K stars, were not an appropriate choice for 
the cold carbon stars. Therefore, the distinction bet ween an R8, the redder R -type stars, 



from an Na, the most blue N-type, was not obvious. iKeenan fc Morgan! (1194 ll ) decided to 
rearrange all the carbon-rich stars in another classification scheme that presented a single 
temperature sequence. The numerical temperature type was, at this time, determined 
based on the atomic and molecular structures, considered more susceptible to temperature 
variations. They also added a C2-index, based on the strength of the C2 bands. Many 
attempts were made to improve this MK system, while keeping its basic structure. 



Yamashita (1972 



19751 ) with the C-Classification System listed in his tables a very large 



-4- 



number of stars with the two main parameters from the MK System, and additional 
intensity indices of several other atomic and molecular characteristics from carbon stars. 
The result of this attempt, was a detailed notation for carbon-rich objects, but it was not 
very practical or compact. 



Later, iKeenanl (119931 ) pointed out several reasons for replacing this old C-Classification 
System, since much importance was given to the Na D-line in determining the numerical 
indices of temperature. The assumption that the Na D-line could be a good tracer of 
temperature in this case, however, was revealed to be quite unfortunate, specially in the case 
of N-type stars, as they have an enormous molecular opacity in the same spectral region 
due to the (7,2) CN band. Furthermore, the N-type and the R-type stars, in fact, describe 
two different populations that should not be classified under the same spectral class. The 
new classification system proposed by Keenan93, the MK Revised System, re-established 
the spectral subclasses for the carbon rich objects and the temperature indices based on 
infra red intensities. Ad ditionally, this Revised MK System listed four abundance indices 



from lYamashital (119721 ): the intensity of the C2 band, the isotopic carbon ratio, the SiC2 



band and the CH band strength. These indices and system are the most w i dely accepted 



and used nowadays in 1 


;he study of carbon stars ( 
)• 


Hatzidimitriou et al. 


2003: 


Morgan et al. 


2004; 


Deroo et al. 


2007 





Concerning the circumstellar analysis, 



Sloan et al. 



(I1998I ) presented the first 



classification based on an infrared study of the spectra of carbon and oxygen-rich AGB 
stars. Qualitative classes for the circumstellar material were settled through two indices: 
one chemical and other a strength indicator of the characteristic structure of the emission 
spectrum. 

In the present paper we present a new scheme of classification for C-rich AGB stars 
(hereafter NSCC) that includes a notation conceived through a wide analysis of the 
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main spectroscopic features of these carbon-rich objects with circumstellar material in its 
surroundings. The main difference between this scheme and the others available is that, 
instead of using a single region of the spectrum, we suggest a classification system based on 
a large wavelength range, from the blue-optical to the mid-infrared. Thus, this notation 
draws a more complete scenario of each star. The parameters used in this extended 
classification system were devised from the analysis of a sample of 51 observed AGB carbon 
stars. 

Section [2] describes the optical observations and other infrared data used and section 
E] and H] details how indices and parameters of classification were obtained. While, section 
[3] explains the optical indices, section @] contains the infrared ones. Section [5] shows an 
analysis of the classification indices, a discussion of possible evolutionary sequences for the 
AGB carbon stars and also an atlas of optical spectra and infrared radiative transfer models 
of the stars studied in this work. The overall results are discussed in section |6j 



Observations and Other Data 



Most of the stars observed were selected because either they have a dubious or 
incomplete classification or doesn't even have one defined. Our intent was to study these 
targets in a large spectral range and, therefore, st ars with data on the Atlas of Low 



Resolution Infrared Astronomical Satellite Spectra (jOlnon et al 



1986|, her eafter LRS-IR AS) 



or on the Short- Wave Spectrometer of the Infrared Space Observatory (IKessler et al. 



1992 



hereafter SWS-ISO), were preferred. In few cases, the selected target lacked any 
available infrared data. We aim to observe those in this spectral region to complete their 
classification, in the next years. 



Another 6 targets, W CMa, RY Mon, V Hya, TW Oph, RY Hya and BE CMa, were 
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used as primary standards in optical range and for the temperature index. Th ese enable a 
calibration of this new scheme here presented with the old ones selected from 



( 11972 



19751 ) 



Barnbaum et al. 



(11999 1 and 



Yamashita 



Keenanl (119931 ). Although some other targets 



from our sample have been studied by these authors, the 6 standards selected either figure 
simultaneously in more than one of the three catalogues or were well studied by others. 

The sample comprises 51 stars which can be seen from the South Hemisphere, that are 
bright enough and not too cold to be observed with an appropriate signal-to-noise ratio. 
Most objects of the sample, 42 targets, were optically observed at the 1.52m telescope at 
the European Southern Observatory, La Silla in several observational runs between 1996 
and 1998. It were used the Boiler and Chivens spectrograph facility with a Cassegrain 
f/14.9 focus. Some of the gratings used at that time are no longer operational. The #11 
grating was centered on 6000A with a dispersion of 66A/mm. The #23 had the central 
wavelength set to 5600A with 126A/mm dispersion. Grating #26 was centered on 5800A 
and 5900 A with 66A/ mm dispersion. Finally, the #32 grating had the central wavelength 
at 5000A and 23A/mm dispersion. The higher resolution was achieved with the #32 with 
0.48A/pix, followed by the #11 and #26 with a 0.99A/pix and the #23, with 1.89A/pix. 
The spectral ranges and dispersions were selected in order to obtain as much details as 
possible in low resolution spectra. Several objects were observed with at least two different 
gratings and, of those, some were observed in four. 

Eight more stars were observed with the 1.60m telescope and the Coude spectrograph 
at the Observatorio do Pico dos Dias, Itajuba - Brazil, during 2006 and 2007. We used 
the CCD #098 which provides a low resolution spectra of 0.25A/mm. To cover the same 
spectral region of the ESO spectra, two regions were observed: A 4300A - A 5300A and 
A 5700A - A 6300A. 



All optical spectroscopic data were post-processed with the Image Reduction and 
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Analysis Facility (IRAF), including wavelength and absolute flux calibration. For the 
wavelength calibration, He-Ar and Th-Ar lamps images were used as reference, which 
were also observed at La Silla and Pico dos Dias, respectively. Several spectrophotometric 
standards, such as HR1541, HR3454, HR4468 and HR5501, were observed to allow flux 
calibration. Both calibration tasks are critical for all the measurements in this work, and 
these data reduction enabled a detailed analysis on the optical spectral features of each star 
of the sample. 



The infrared spectra used were taken from ISWS-ISO 



and 



1 

LRS-IRASl databases. These 



two catalogues together give a spectral range from 2.4 to 45 /mi. The former provides 
a better resolution (circa 1500 - 2000) than the latter (circa 20 - 60), as well as a wider 



spectral range. Unfortunately, only two object s from our 



For the other stars we had to rely only on the 



sample have an 



SWS-ISO 



LRS-IRASl spectra. 



spectra. 



3. Optical Classification 

The notation ado pted in this w ork for the optical properties is the one described on the 



MK Revised System ( iKeenan 



19931 ) with some slight improvements. The enhancement at 
the notation is simple and self-explanatory, and its four indices allow a better sketching of 
the scenario of an AGB carbon star photosphere. The luminosity classes were not applied, 
as they are not well defined for most of the objects. As usual, an uncertainty character (:) 
was added, and it suits indices that were calculated with a poor quality spectrum. 

Most optical parameters are measured relative to the continua, but unlike G- and 
K-type stars, to establish a continuum for an AGB carbon-rich star is not a trivial task, due 
to the high opacities on its optical spectrum. Therefore, local continua or pseudo-continua 
were used instead, defined by regions where C2 e CN absorption bands are weak. As a 
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result, some uncertainties to the parameters are added, that are not directly measurable. 
Nevertheless, all the parameters used to calculate our indices were obtained by averaging 
measurements from multiple observations of each star, when it was possible. At least two - 
and in several cases four - different spectra were used in this average, decreasing the effects 
of uncertainties in the measurements. 

Only the three first indices, C2-, j- and MS-index, had their parameters measured 
relative to the continua, and for those, two pseudo-continua were drawn: one for the first 
two indices and another for the later. The first pseudo-continuum was defined as a legendre 
function fit between maxima at A 5722A, A 6202A and A 6620A, as can be seen in Figured] 
(left). These are points where the absor ption band of C2 and CN is weakest and have been 



suggested for AGB carbon-rich stars by 



Westerlund et al 



(09911) and 



Hatzidimitriou et al. 



( 20031 ). For some stars, the maximum at A 6180A was greater than the one at A 6202A, 
so the pseudo-continuum was defined without that point. The second, Figured] (right), 
was a linear local continuum defined between maxima at A 4962A and A 5030A. These 
maxima were suggested by 



Sarre et al. 



(120001 ) because no overlap from metallic lines and 



other bands compromise the fluxes at the points. All this was done after applying a slight 
smoothing to the spectra using a 3-pixel box and then, the spectrum was normalized 
related to the pseudo-continuum settled. The strengths were calculated with respect to the 
pseudo-continuum and also to a local continuum defined by a nearby maximum. 



3.1. The Comdex 



The more essential abundance index is the C2-index, which is an index of the carbon 
excess over the oxygen in th e carbon-rich photospheres. This index was originally defined 



in the Revised MK System ( IKeenan 



19931 ) as dependent in strength of the Swan Bands A 



5165A and A 5635A. However, these bands are just suitable for C-R and C-H, since in these 
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types of stars the C2 bands are strong; but for the N-type and most of the J-type stars, 
that have a compromised flux at this spectral range, it is difficult to obtain the C2-index 
through analysis of these two specific bands. Fortunately, many Swan bands are widely 
present at the AGB carbon star spectra, and several may be used for the same purpose. 

Thus, we decided to adopt the well correlated parameters suggested by 



Hatzidimitriou et al. 



(120031 ) in their study of a LMC carbon sample: the sum of the 
strengths of the normal and isotopic C2 band, associated with A 6192A and A 6168A 
respectively, shown in Figure [2] (left); and the equivalent width of a complete absorption 
in the range of A 5722A - A 6202A, shown in Fi gure[2] (right). The strength of the bands, 
Da 6168 an d D A 6192 , were measured from the minima, at the wavelength described, with 
respect to the pseudo-continuum for normalized spectra. The complete absorption by 
C-rich molecules associated with the equivalent width, Wa 5722 - a 6202, were measured in 
the edges at the maxima described. Both parameters, strength and equivalent width of 
the bands, are measurements of the same physical property, i.e. the carbon abundance of 
the photosphere. That means, these are expected to have a positive correlation, since an 
enhance on the carbon abundance will reflect in a more strong strength of the bands as well 
broader equivalent width. 

It is also important to notice that the isotopic band is significant only in stars with 13 C 
excess, and for those it has to be considered. For the ordinary carbon stars, the spectrum 
between A 6159A and A 6175A shows a broad C 2 absorption feature with a minimum at A 
6162 - A 6164A and a secondary minimum at A 6168A. On the other hand, the stars with 
13 C excess show a single minimum at A 6168A. An extra difficulty in identifying 13 C 12 C A 
6168A, when it is weak, is the proximity of Cal A 6162A and several CN and 12 C 12 C bands. 

The former parameter, D A 6192 + D A 6168 , is obtained in a shorter spectral range, and 
therefore it is a more reliable measurement, once the local continuum is better defined, and 
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for this reason it was designated as the primary parameter, p\. Thus, the equivalent width 
of the large range, A 5722A to A 6202A, was designated secondary parameter, p 2 . If, 
for a given star, only pi can be measured, it is usually reliable enough to be used alone, in 
place of the average indicator described bellow. 

This hierarchy of trust were used, then, to re-scale p 2 , by putting it at the same scale of 
the primary parameter. Some stars studied have multiple spectra and, for those an average 
of each parameters was calculated and it was also possible to assign an error for p\ and 
P2, caused by the use of different telescope and grating combinations for the observations. 
Then, after measuring the both parameters, a linear least-squares fit, p 2 = a ■ Pi + b, was 



calcu 



Hatzidimitriou et al. 



Keenan 



ated, wielding a = 242 and b = 38. The fine correlation found by 
(2003) can also be seen for our Galactic C stars, as were suggested previously by 
( 11993I ). and so the linear fit were adopted. The coefficients a and b were used to apply a 
linear transformation to the secondary parameter, in order to get a new one, p' 2 , in the 
same scale of p\. 



/ P2 

P 2 = — 



(1) 



The next step was the evaluation of the average indicator: 

- P2 + Pi 
P = — o — • 



(2) 



This indicator is used to estimate the ranges of the levels of the carbon abundance index. 
This simple average, between the parameters previously set to the same scale, allows a 
further reduction in the uncertainties of the measurements. 



In 



Keenanl (119931 ). the C2-index has a increasing scale of eight levels, varying from 1 to 



8, which we adopted. And the range values given by the average indicator for each index 



l evel were se t by 



(11972 



comparing y of t h e 6 st andards objects to the C2-index set in 



1975|) and 



Barnbaum et al 



Yamashita 



(119961 ). None of these schemes have given the boundary 
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values of their domains, nevertheless, we were able to make our scale agrees closely with 
the index of excess of carbon over oxygen given by Yamashita's tables. The ranges of each 
level of the C2-index , which have an equally interspace, is given by the average indicators 
described on Table [H 

The C^-index can be calculated for any carbon star, by obtaining pi and P2 from its 
spectra, then using eqs. [1] and [2] to obtain p and consulting the Table [H The results 
of the C2-index obtained for our entire sample are discussed in Section [51 Therefore, the 
C2-index, in this way described, corresponds to a direct measure of the carbon intensity 
observed on the carbon-rich star, and not one relative to a standard of same temperature, 
as used in the past. 

As multiple spectra for some stars were taken, including different combinations of 
telescope and gratings, it is possible to expose how consistent is the C 2 -index of this new 
scheme for a same star observed with different resolutions. We manage to have different 
combinations for 20 targets. The parameters p\ and P2 were obtained for each star spectra, 
and then, a deviation from the average value for each parameter of each target were 
calculated. In order to get a global behavior, the parameter's deviations were averaged 
through the whole sample and it was found an error of 5 P1 = 0.06 and 5 P2 = 13.80 for the 
primary and secondary parameters, respectively, which means an error of 5p = 0.04 for the 
average indicator. These errors, compared to the magnitude of p\ and P2 obtained and p 
calculated, indicate a good reliability in the method and in the bands adopted. In fact, the 
strongest C2 bands of the Swan System can be estimated even on prismatic spectrograms of 



very low resolution ([Blanco et al. 



1978 



Westerlund et al. 



1978 



Sanduleak fc Philip 



1913). 



On the other hand, a discussion about how consistent our C2-index is with previous schemes 
is placed in section El 
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3.2. The j-index 



Let us discuss now the isotopic carbon abundance. It deserves a specific index, because 
of its historical importance and it is extremely releva nt in stars with low 12 C: 13 C ratio. 
The expression J-Star, used for the very first time by iBouigud ( 119541 ). refers to stars that 



have an unusual strong isotopic ca rbon bands . The exact evolutionary st age and nature of 



this C-rich stars, is still unclear (Abia et al. 



20031 ). since 



Utsumi 



(1985) found no excess 



of s-process's elements in J-type stars, as it would be expected for any ordinary AGB 
carbon star. An inspection in R-type stars reveals that th ey also do not present an excess 



of these elements together with an overabundance of 13 C (iDominy 



1984|). Despite both 



C-R and C-J may seem to have their carbon -rich photosphere orig inated with the same 



mechanism, with an anomalous helium-flash ([Dearborn et al. 



20061 ). no connection between 



both types of stars was found yet. So it is not certain whether C-J evolved and cooled 
from the C-R. The absence of s-process's element enrichment does not fit these stars in the 
AGB scenario that supposedly already crossed the thermal pulses, discarding th e possibility 



of a 13 C enrichment due to the third dredged-up (IPilachowski Sneden 



1988|). The Tc 



detection together with the s-process's element absence in two C-J Peculiar Stars made 
them even more riddling, as the Tc presence points dredge-up occurrence. Some attempts 
have be en made to explain this apparent incoherence through a cool bottom processing 



theory (IWasserburg et al. 



19951 ). However no conclusion about the nature of the J-Stars is 
presented, these objects are represented in the NSCC by a C-J notation for the spectral 
sub-class index. We intend to stress their unusual nature and to distinguish them from 
stars with normal 13 C abundance. 



Hatzidimitriou et al. 



(120031 ) found a strong correlation between the ratio of the strength 
of the isotopic and normal C2 band, associated with A 6168A and A 6192A respectively, 
shown in Figure [2] (left); and the ratio of the isotopic and normal equivalent width of the 



13 



CN bands whose band-heads are located at 13 C 14 N A 6260A and 12 C 14 N A 6206A, shown in 

;he isotopic an d 



Keenanl (119931 ). in 



Figure [31 for LMC stars. These strength ratios and equivalent width of 
normal bands are the ones suggested in the study of Galactic C-Stars by 
a previous work. In our scheme of classification these parameters were also used to evaluate 
the isotopic carbon abundance, namely the j-index. As done for the C2-index, the strength 
of the bands, D A 6168 and D A 6192 , were measured from the minima of the band with respect 
to the pseudo-continuum at the normalized spectra, and the equivalent width, Wa 6206 an d 
Wa 6260) were measured at the maxima of the edges of each waveband. Regarding the 
isotopic and normal CN bands for C-J stars, it can be seen a clear splitting of the waveband 
in two halves A 6205 - A 6250A and A 6255 - A 6290A; while for normal N stars this splitting 
are not so obvious. 

All measurements for the j-index are taken in short wavelength ranges, being equally 
reliable. The ratio of the strengths of the C 2 bands, D A 6168 / D A 6192 , was designated as 
parameter qi, and the ratio of the CN bands, Wa 6260 / Wa 6206, as parameter q 2 . For the 
stars of the sample that have multiple spectra were possible to assign an average for each 
parameter and an error associated with the effects in qi and g 2 caused by the use of different 
telescope and grating combinations on the observations. 

Both parameters are measurements of the 13 C excess, thus it is expected that they 
have a positive correlation, since an enhance on 13 C will reflect in both measurements. 
Nevertheless, a fine correlation can only be seen for C-J stars. For this reason, only 
likely J-stars could be used to calibrate the av erage indicator formula. We selected those 



candidates by applying the definition given by iGordonl (119681 ) for J-stars, as the ones that 
have the isotopic band 13 C 12 C A 6168A with at least half the strength of the normal band 
12 C 12 C A 6122A. 

The qi and g 2 obtained for the targets selected through Gordon's definition seem to 
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have a linear correlation as were suggested by iKeenanl (119931 ). Then, a linear least-squares 
fit, q 2 — c ■ qi + d, was calculated for this small sub-sample, wielding c = 3 and d = —1. The 
coefficients c and d were used again to apply a linear transformation on q 2 to acquire two 
parameters at the same scale, q' 2 and q%. Thus, a simple average between these two gave the 
average indicator, q. 

From the targets selected through Gordon's definition only 3 stars of our sample 
(BE CMa, BM Gem and SU Sco) already had a j-index available, then, they were defined 
as the standard objects. Following the MK Revised System, the j-index were divided 
into eight l evels, varyi n g from to 7, which also agree with the isotopic carbon strength 



index from 



Yamashital ( 1972 



19751 ). As presented in Section I3TTI the ranges of the average 



in dicator for each le vel w ere settled by comparing the q of the standards to the j-index data 



m 



Yamashital (119721 ) and 



Barnbaum et al. 



(119961 ) . Again, the boundary domains of each 



catalogue were not described but we were able to make our scale agrees closely with the 
index of 13 C excess in Yamashita's tables. The ranges of each level of the j-index , which 
also have an equally interspace, is given by the average indicators described on Table El 



Usually, C-N stars have j-indices between 1.5 and 3 and C-R between 2.5 and 4 (IKeenan 



19931 ). Thus, in the NSCC, only targets with j ^ 4 are surely considered a J-star and have 
an explicit j-index. In all other cases, the value of the j-index is omitted. The j-index can 
be assigned for any carbon star, by obtaining qi and q 2 from its spectra and consulting the 
Table [2j The results of the j-index obtained for our entire sample are discussed in Section |5j 

Through the multiple spectra, with different combinations of telescope and gratings, 
taken for 20 targets, it were possible to expose how consistent is the j-index of the NSCC 
for a same star observed with different resolutions. First, a deviation from the parameters 
qi and q 2 were calculated as done in the previous Section. Then the global behavior, an 
average of the parameter's deviations through the hole sample, was evaluated and it was 
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found an error of 5 qi = 0.03 and 5 q2 = 0.09, which means an error of 5q = 0.02 for the 
average indicator. These errors, compared to the magnitude of q\ and q<i obtained and q 
calculated, also indicate a good reliability in the method and in the bands adopted. 



3.3. The MS-index 



Merrill! (119261 ) and ISanfordl (119261 ) observed for the first time a molecular absorption, 
in the violet-blue region, due to a compound associated with the SiC 2 molecule. Some 
SiC2 absorptions (the Merrill- Sanf or d bands, hereafter MS-bands) can be found in the 
A 4300A - A 5000A spectral range. The MS-bands are not present in all C-rich stars, but 
when present, they are very strong (see Figure HI) and, therefore, a representative index of 
the stre n gth o f the SiC2 molecule in their photospheres can be obtained as suggested by 
Keenanl ( 119931 ). The C-J stars offer favorable conditions to their study, since they usually 
present MS-band structures, and their flux at the violet-blue range are not as compromised 
as one of a C-N star. 



The MS-index is the last optical index added to the NSCC notation for the description 
of the photosphere of the carbon-rich stars. Unfortunately, not all stars in our sample were 
observed in the spectral region populated by SiC2 structures, thus, it is possible that some 
targets, which have a strong MS band, do not have a MS-index in the final classification. 
As it is a complementary index, we considere d that there i s no need to subdivide it into 



too many levels as it has been usually done (Keena 



in 



19931 ). Our scheme of carbon star 



classification presents the s trength of t 
medium and strong (as in 



Sarre et al 



re Si C2 index divided into only tree levels: weak, 



20001 ) . The parameter employed to obtain the 



MS-indices were the strength related to a local continuum, d, of a SiC 2 structure with 
head-band located at A 4977A. 
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The rang e s of e ach level of the MS-index in the NSCC are those described by 



Morgan et al 



(120041 ) . that also used the same band for d, and vary from 1 to 3 in a 
increasing scale. The NSCC has an additional MSO level, for very weak or absent MS bands, 
which can be omitted from the final notation. 

Using different combinations of telescope and gratings, taken for 16 targets, it were 
possible to expose how consistent is the MS-index of the NSCC for a same star observed 
with different resolutions. A deviation from the parameters d could be calculated for these 
stars, and then, the global behavior error, an average of the parameter's deviations through 
the hole sample, was evaluated in 5d = 0.03. 



3.4. Spectral Subclasses 

Due to the high opacity in the optical spectra of the cool carbon stars, several criteria 
have been suggested to define each spectral subclass for carbon stars, but each single one 
did not allow by itself a consistent classification. In this work, we made an inspection 
through all visible spectrum, by combining the use of some of these criteria. When applied 
together in a simple yes/no logical analysis for each characteristic feature, they complement 
each other, giving a more conclusive result. 

The first criterion concerns the extension of the spectra to the blue-violet range. 
Although a high opacity in wavelengths lower than A 4977A points to a carbon C-N star, it 
is not a sufficient condition, as cool C-R and C-H stars may also exhibit this effect. 

The presence of structures related to s-process elements are easy to identify in C-N and 
C-H but not in C-R stars, which were not enriched with these elements. A good method 
to qualitatively detect the presence of these structure lines is: when the intensity of the 
Ball lines A 4554A and A 4934A is stronger than the closer CN bands, for instance, the 
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band-head at A 4576A; an d when the Sri A 4607A has about the same intensity of the band 



this might be a C-N star (IBarnbaum et al 



1999 1. In FigureE] (left), the RV Cen star has 
the Ball and Sri lines stronger as the CN band, but the same does not happens for C* 1130. 

Finally, a third inspection concerns the distinction between C-H and C-N stars, as both 
present strong Ball lines with respect to the CN band. It is possible to distinguish these 
type of carbon stars by presence of the P-branch, which band-head is located at A 4352A. 
When P-branch is evident, this is typical of C-H's and not expected for C-N's, as it can be 
seen in Figure [5] (right). This branch is also present in most of the C-R stars. A summary 
of the combined use of these criteria can be seen in Table [3J 

All these inspections (the extension to the blue- violet spectral range, the presence of 
s-process element's lines or the presence of the P-branch, A 4352A) are just a qualitative 
analysis based on a visual inspection. On the other hand, such a study can be applied 
even to prismatic spectrograms of very low resolution, as no high precision is required 
to establish the spectral subclass through this method. Anyway, these three inspections 
together can give an unambiguous way to assign the spectral subclass of each carbon star 
studied. The spectral subclass goes at the beginning of the NSCC notation. 



4. Circumstellar Envelope Classification 

During the Asymptotic Giant Branch phase, stars suffer intense mass loss which leads 
to the formation of circumstellar dust envelopes. In order to describe a complete scenario, 
three more indices were added to NSCC: two concerning the circumstellar environment 
itself and one related to the temperature of the star. This infrared classification, which 
we are introducing in this work, was obtained through a Monte Carlo numerical algorithm 
which simulates the radiative transfer problem in the circumstellar envelopes. 
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The numerical treatment of radiative trans fer applied to a spherical envelop e, used 



Lorenz-Martins &: Lefevrd (I1994J ). Here 



by us here, was described in a previous work of 
we present it briefly for the sake of completeness. The propagation of stellar and grain 
radiative energy is simulated photon by photon following a Monte Carlo scheme. For each 
interaction between a photon and a grain, a fraction of the energy is stored (absorption) and 
the remaining part is scattered according to the scattering diagram. The stellar radiation 
leads to a first distribution of dust temperature and the thermal radiation from grains is 
simulated, giving after several iterations the equilibrium temperature. 

When two kinds of grain are present, the path length of photons is defined along any 
direction by the total opacity 

r ext (X) = Ti(A)+T 2 (A) 

where, 

7i(A)= / Ni(r)rra^Q exU (X)dr 
Jo 

in which a« is the radius of the grain i, Q e xu{^) its extinction efficiency and iVj(r) its 
number density at the distance r from the center of the star. Efficiencies and scattering 
diagrams are computed at each wavelength using the Mie theory. When the position of the 
interactive grain is defined its nature is determined: the probability of having a grain % is 

Pi{\) = a*Q ext Mr)/Y. a2 ^ N ^- 

j 

The energy absorbed is then calculated with the relevant Q a bs and the new direction 
of scattering is generated according to the corresponding scattering diagram. The energy 
is stored separately for each kind of grain and then two sequences of temperature are 
obtained. The origin of the emergent radiation at any wavelength is easily identified (direct 
stellar radiation, scattered light, emission from each kind of grain). 

The mixture is characterized by the ratio Ni(r)/ N 2 (r). Here, this ratio will be 
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independent of the position point. When one specie is much less abundant than the other, 
a large number of events have to be simulated to get the good statistics. The total number 
of grains is defined by the extinction opacity calculated along the radial direction between 
the inner radius of the shell Ri and its outer radius R Q at A = 1/im. The number density is 
assumed to vary as r 2 , corresponding to an expansion at constant velocity. Therefore, the 
following physical quantities are required: 

• the effective temperature of the central star T e ff] 

• the inner and outer radii of the shell Ri and R a . We shall assume here that both kinds 
of grain exist within the same limits. This point will be discussed later; 

• the grain radii and refractive indices at all wavelengths. (A grid of 30 wavelengths is 
used) ; 

• the ratio N x / N 2 ; 

• the extinction opacity at A = 1/xm. 



The computation gives the spectral distribution of the total flux and of its different 
components (direct, scattered, emitted), the temperature law for each kind of grain. Best 
fits ive the effective temperature of the central star, SiC/A.C. ratio, r and grain size. 

The absorption and scattering efficiencies, as well as the albedo for the grains were 
calculated by us using the Mie theory and the optical constants tabulated in the literature. 
For the amorphous carbo n (hereafter A.C.) we have used the optical constants published by 



Rouleau &: Martini (1l99fl ) since they were calculated for a wide range in energy (4.127 ■ 10 s 



eV - 3.500 • 10 3 eV) satisfying the Kramers- Kronig rela tions. For SiC g rains we used the 



optical constants of silicon carbide (SiC) d etermined by 



those obtained by 



J5) 



Pegouriel (119881 ) and for silicate 



David h Pegouriel (119951 ) . These last were used only to model silicate 
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carbon stars. Finally, the infrared notati on adopted is based on paramet ers and ranges of 



each index level proposed previously by 



Lorenz-Martins &: Lefevrd (119941 ) 



Although empirical measures of the size of the emission feature could be applied with 
some adjusts in the methodology here presented, it is not advised. A nyone intere sted in the 



19961 ). and the 



NSCC can use a public radiative transfer code, for instance DUST (jlvezic 
quantities needed in the scheme can be obtained. 

The circumstellar envelope classification completes the NSCC, and can be filled using 
any infrared spectra that provides wavelength ranges around the dust features, i.e. 11.3 
fim for normal carbon stars and 9.8 /im for the peculiar ones. There are about 900 carbon 
stars observed by IRAS satellite and almost 50 by ISO satellite. Furthermore, it is possible 
to obtain new fresh IR data with Spitzer satellite and, at ground, using Gemini telescopes 
with TReCS and Michelle spectrometers which provide wavelength ranges required for the 
NSCC indices. 



4.1. The r-index 

The r-index or opacity index is an estimate of the circumstellar envelope's optical 
depth, and can provide information about the evolutionary stages of the carbon-rich stars, 
as discussed in section disc. Through the opacity index it is possible to establish how 
optically thick is the envelope, allowing us to draw a better scenario for the star-envelope 
system. 

The parameter used to define the r-index is one of the outputs of the radiative transfer 
model fit, the r parameter. We divided this index into four l evels, varying from to 3, 



and th e ranges of each level agree with the ones described by 



Lorenz-Martins fc Lefevre 



(jl994j), as can be seen in Table HI The rO level corresponds to an extremely low 
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opacity value, nevertheless , both rO and rl correspond to the opacity level I, defined by 



Lorenz-Martins fc Lefevre 



The deviation of the optical depth found through the radiative transfer code applied 
to a spherical envelope was S T = 0.03. T his is a zero-point error for all the circumstellar 



envelope's optical depth obtained with the 



Lorenz-Martins fc Lefevrd code. When compared 



5 T to the ranges of each level of the r-index, we verify the good reliability in the method 
adopted. 



4.2. The SiC-index 

The SiC-index is a relative abundance of silicon carbide to A.C. (hereafter SiC/A.C.) 
for the grains in the circumstellar envelope. This ratio abundance is based on the spectral 
feature at 11.3 /im, which appears in most of the carbon stars. Modeling their circumstellar 
envelope by taking into account only the amorphous carbon dust grains is a good first 
approximation, since it corresponds to the main emission of the envelopes. Nevertheless, it 
does not reproduces the emission feature at 11.3 /im, associated with the silicon carbide. 

Therefore, the dust compounds applied to obtain the SiC-index were the SiC/A.C. 
ratio abundance, which is also an output parameter of the radiative transfer model. A 
model fit to the spectral feature at 11.3 /im can be seen in Figure [6] (a), in dashed line. 
The open circles in the figures denote the photometric data taken fro m SIMBA D database, 



SWS-ISOl data. The 



ranging from optical to infrared wavelengths, while the full line is the 
deviation of the SiC/A.C. ratio abundance found, for all models calculated, through the 
radiative transfer code was Ssic/A.c. — 0.02, which also represents a good trustworthiness in 
the method. 



A simplified SiC index was presented by 



Sloan et al. 



(119981 ) . They have developed an 
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infrared spectral classification for 96 carbon star based on iLRS-IRASl spectra. Basically 
they extracted the blackbody contribution using a 2400 K Plank function for their sample. 
One problem is that the residual spectra include not only the dust features but also the 
photospheric absorption bands. In addition, they used a flux ratio to obtain their classes. 
Even so this is an interesting way to classified dust envelopes, we believe that both r- 
and SiC- (or Silicate) indices are better obtained using a full modeling including radiative 
transfer. In this case, we can estimate the contribution of any kind of dust present in the 
circumstellar envelope. 

We set four levels for the SiC-index, v a rying from 1 to 4, and the domains of each 



level are listed in the Table [51 



Sloan et al. 



(119981 ) described three indices related to the 



intensity of the emission feature and others related to the surrounding continuum and to 
the unusual profiles found by them. When the levels o f the SiC-index are compared with 



the description of the first three indices of 



Sloan et al 



, we have, roughly, that: our SiCl 



and SiC2 correspond to their SiC class, our SiC3 to their SiC+ class and our SiC 4 to 
their SiC + + class. 



The dust species formed in the circumstellar envelope of car bon stars shou 



reflect the chemica 



compo sition of their photospheres. H owever. iLittle-Marenin 



LRS-IRAS 



d, at fi rst, 
(I19861 ) and 



catalogue, discovered 



Willems &: de JongJ (119861 ). independently, inspecting the 
the existence of a sample of carbon-rich stars which exhibit oxygen-rich dust envelopes. For 
these peculiar carbon stars, the carbon-rich nature is revealed in the visible spectral region 
while the oxygen- rich o ne is exposed in the mid-in frared, where amorphous and crystalline 



silicates can be detect. 



Willems fc de Jongi (119861 ) showed that these stars are, in fact, 



J-type carbon stars. In our scheme of classification, these silicate carbon stars received 
a J Pec notation instead of a SiC-index. Figure [6] (b) presents the best fit model to the 
spectral feature at 9.8 /j,m associated with silicates, typical of a silicate carbon star. 



-23- 



4.3. Temperature Index 



It was only with iKeenanl (119931 ) that a system of classification of carbon stars 
incorporated a numerical index of temperature based on infrared fluxes. Those indices were 
more reliable than the atomic lines commonly used to estimate the temperature for stars of 
the class M and hotter. 

In the NSCC, the temperature index is also devised from a parameter based on infrared 
fluxes by way of the radiative transfer simulation. Once the modeled SED which best fits 
the infrared data is set, the effective temperature of the central star (T e ff) is one of the 
outputs of the simulation. And that is the parameter that should be used to establish 
the temperature index levels of each target. The temperature index is placed, in the final 
notation, just after the spectral subclasses, as usual. 

By means of this radiative transfer code applied to a spherical envelope we were 
able to settle an error nearly 100K for the effective temperature of the central star 
calc ulated. This is a zero-poin t error that should be considered for all data modeled with 



the 



Lorenz-Martins fc Lefevre code. Above this error value, the calculated models differs 



meaningfully of each other. 

The range values of the effective temperature for each index level w ere set by compar ing 



this output pa r amete r of the 6 standards objects to the indices set in 



1975|); 



Keenan 



(11993f) and 



Barnbaum et al 



Yamashital ( 1972 



(119961 ). These boundaries can be seen in Table 



[6j Even knowing that some of these previous schemes employed parameters others than 
based on infrared fluxes to set their temperature scales, we made an effort to agree our 
indices levels with these old ones, by just reformulating the range values of each level of the 
temperature index. 
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5. Sample Classification 

As a first application of the NSCC, we calculated all the parameters presented in 
sections [3] and H] for all stars in our sample. Thus, up to seven indices of classification were 
established for each target, including for the 6 standards used in the development of the 
methodology, that can be seen in Table [HJ 

5.1. Optical Parameters 

Through the spectral subclasses we identified not only C-R and C-N stars in the 
sample but also four C-H candidates, which should be further investigated, to a final 
conclusion. Among those, W CMa and V Hya have currently a C-N classification, and as it 
was not possible to analyze other reliable CH bands, they received a C-N: index in the final 
notation, denoting an uncertainty. Some stars previously classified as C-N, were identified 
as actually belonging to the C-R spectral subclass: DH Gem and UW Sgr were reclassified 
from C-NO to C-R and NP Pup from C-N4.5 to C-R. 

The spread of our sample in regards to the parameters p\ and p' 2 for the C2-index and 
the established ranges of the average indicator, p, can be seen in Figure [7] (left). The full 
line illustrates when primary and secondary parameters match, while dashed lines mark the 
edges of each level range. It's also possible to see, in this figure marked with open circles, 
the standard objects, which were used to obtain the ranges of each C2-index level. The 
stars with the higher C2-index, namely C26 and C27, were named in Figure [7J (left) and T 
Mus has the higher estimate of carbon abundance. On the other hand, HD 113801 and GP 
CMa have the lower C2-index calculated, C2-2. 

When the j-index w as assign ed to all stars, it revealed 14 targets as J-stars, three more 
I I 

than if selected only by iGordonr s definition: V971 Cen, CL Mon and C 1130. Figure [TJ 
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(right) displays the edges of each level of the index in dashed lines and plots the j-index 
obtained for the whole sample. The correlation between the parameters q\ and q' 2 exists 
only for J-like stars, which means that, this behavior is an evidence of stars with 13 C excess 
only. It can be seen that, indeed, it do not holds for stars with j ^ 4. As done before, the 
objects used to calibrate the ranges of each level are marked with open circles. The stars 
with the higher j-index, namely j7, are named in Figured (right) and, again, the isotopic 
carbon abu n dance estimate for T Mus is the higher one. We confirmed the classification of 



Chen et al. 



(120071 ) for S C en, V971 Cen and C* 2 208 as C-J stars, in cont rary to the C -R 



Barnbaum et al. 



(|l996|) for the first two and bv lChanl (1199.1 ) for 



classification presented by 
the later, in the past. 

The average parameter d were used to obtain the MS-index for the 32 stars of our 
sample. Comparing the MS-index with C2-index it is possible to see that the stars with 
the higher value of the SiC2 abundance, index MS3, seems to be also those with the higher 
values of C 2 -index (C 2 6 and C 2 7): S Cen, TZ Car, W Pic C* 2208, FO Ser and T Mus, for 
instance. 

The optical spectra of the stars classified in this work can be seen in Figures [HJ to [T5l 
The most significant features for the classification of these objects under the NSCC are 
marked on the spectra whenever present. The spectra are displayed in a decreasing scale of 
temperature separated by their spectral subclasses: C-N stars are shown in Figures [HI to [T2| 
C-R in Figure [13] and C-J in Figures [TH and [15j C-H stars are displayed together with the 
C-N due to their spectral subclass uncertainties. 
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5.2. Circumstellar Envelope Parameters 

We also have determined, for the first time, circumstellar envelope opacity index 
for 37 stars and SiC/A.C. ratio abundance index for 26 stars. C* 1003 and HV CMa 
obtained the higher values on r-index, suggesting an optically thick envelope usually 
associated with tip-AGB stars. It is possible to assume, that these stars already ejected 
a considerable amount of matter, which condensed and accumulated, creating a denser 
envelope. Nevertheless, some stars have an extremely low opacity value, rO, which 
represents the opposite scenario (e.g. NP Pup and BE CMa). These stars have an optically 
thin circumstellar envelope which may correspond to early stages of carbon AGB. In regards 
to the J-type carbon stars, we notice that the star T Mus not only have the higher optical 
indices, but also the higher opacity index, t2. It is possible, than, T Mus may be the most 
evolved J-type star f rom this sa mple, while, on the other hand, BE CMa could be the least 



evolved one. Not all 



LRS-IRAS 



spectra here presented has a good signal to noise ratio, so, 
in these cases, after the circumstellar envelope classification the object with a poor infrared 
datum received an uncertainty index (:). 

All models calculated for the stars that were classified in the infrared are displayed at 
the Figures [16] and [T71 The stars in Figure [161 ordinary carbon stars, had their radiative 
transfer models calculated to fit the 11.3 fim feature, associated with the ratio SiC/A.C, 
while Figure [T7| shows the models calculated to fit the silicate feature at 9.8 /j,m for the 
silicate carbon stars. Both groups are sorted by an increasing scale of their opacity indices. 

The fit of the model is not always perfect. It can be due to different effects. For 



instance, the data of each target were obtained b y several o 
which means that, the photometric data and the ISWS-ISO 



Dservers in different moments, 
spectra were 



or 



LRS-IRAS 



observed at different time. Knowing that carbon stars are variable stars, it is expected 
that the flux varies with each luminosity phase. Another effect is that, in some cases, the 
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circumstellar envelope may not be spherically symmetric, and therefore, the fits, which 
results from the spherically symmetric models, are not always a good first approximation. 

SiC grains may nucleate closer to the star, at higher temperatures of about 1300 
K to 1500 K, a nd therefore they may be the fir st grains to condense in the carbon-rich 



star envelopes ( 



Lorenz-Martins fc Lefevre 



( IFrenklach &: Feigelson 



19941 ). while A.C. grains condense at 1000K 



19891 ). High values of SiC/A.C. ratio denote, then, early stages of 
the AGB. When more carbon is ejected to the circumstellar envelope, during the thermal 
pulse phase, other compounds begin to nucleate and become more abundant than the SiC, 
for instance amorphous carbon grains. Thus, we can expect that more evolved AGB stars, 
besides being optically thicker, have lowest SiC/C.A. ratios. 

It is also possible to assume that the circumstellar opacity increases monotonically 
with time, since as more evolved the star is, the thicker is its envelope. Differently, a given 
SiC/A.C. ratio decreases gradually during the evolution of an AGB carbon star. Figure [18] 
illustrates both parameter behaviors by the global spread of the circumstellar opacity and 
SiC/A.C. ratio values calculated for the sample. 

We can assume that, by the observed spread, as the envelope opacity grows thicker, the 
SiC/A.C. ratio decreases with a power-law: SiC/A.C. = (3r a . Considering the stars that 
exhibit the greater values of SiC/A.C. at each opacity, i.e. the outer points in Figure [T8l 
we may presume that these stars represent the different evolutionary stages of an SiC/A.C. 
extreme object. Therefore, through an empirical fit to these stars we established a as 
approximately —0.5. 

Varying the (3 coefficient it is possible to obtain several evolutionary sequences 
determined by each initial SiC/A.C. ratio. Clearly, we can not consider r = to compute 
the initial ratio, since at this stage level, the SiC/A.C. value would be indeterminate by the 
lack of matter at the circumstellar envelopes. But, assuming an initial opacity of r = 0.25, 
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which we may consider a stage were it start having a significant density of condensates, it 
is then possible to calculate the (3 value. Thus, the evolutionary sequence expression that 
connects circumstellar opacity with SiC/A.C. ratio is: 

SiC/A.C. = SiC/A 2 °-^ r" 1 / 2 . (3) 

Four examples of evolutionary sequences are plot in Figure [18] with the spread of C-N, 
C-R, and C-J stars of our sample. The evolutionary sequence of a C-J stars is still an 
open question, but it surely does not follow the same evolution of an ordinary AGB carbon 
star. These stars are spectroscopically similar to the C-R stars, which present anomalous 
carbon enrichment durin g the second dredge-up, to such an extent that it is possible that 



19861 ). For fact, as a star go further into the TP- AGB 



the later evolved to C-J (lEvansi 
phase, the 12 C: 13 C ratio increase with the thermal pulses, i.e. decreasing the photospheric 
13 C abundance. We may, then, consider that as a C-J star evolves into the TP-AGB, it 
gradually change their photospheric abundance becoming, at the end, a C-N star. It is 
cautious to point out that not all the C-N stars evolves from a C-J star, most of them 
should f ollow the well accepted sequence of M — > MS — > S — > SC — > C-N (see e.g. 



Herwigi (120051 )). But, if we consider the previous suggestion for the C-J evolution, the 



possible sequence that describes it is C-R — > C-J — > C-N. 

Regarding the circumstellar envelope, we can consider that C-R and C-J star becomes 
carbon-rich earlier than C-N stars, as they began their C-rich mass loss earlier. Therefore, 
their envelope should have SiC/A.C. ratios greater than the C-N stars, since C-R and 
C-J began to nucleate grains before the later ones. Additionally, following our suggested 
evolutionary sequence, it is reasonable to expect that a C-N star has a thicker envelope 
than a C-J, and on the other hand, C-J's envelope should be also thicker than the C-R's. 
Figure [18] illustrates this possible evolutionary sequence C-R — ► C-J — ► C-N as can be 
seem in both the sequences with f3 = 0.02 and (3 = 0.24. Of course, a greater sample is 
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essential to confirm this assumption. 



5.3. Comparison with Other Schemes 



A discussion about how consistent is the NSCC compared with previous schemes is 
important to evaluate its reliability. The seven indices of our new system acquired for 
the sample were qualitatively compared with the most quoted schemes for carbon stars. 
Although most of the indices could be examined, the infrared ones are improvements to our 
scheme and not all of them have been consider by other authors. 

The spectral subclasses obtai ned through th e NSCC methodology co uld be compared 



with two schemes of classification, 



Keenanl (Il993h : 



Barnbaum et al. 



(jl99fih . that had 4 and 



11 stars of our sample alr eady published, r e spect ively. The C2-index were examined in 



regards to three schem es: 



sample, 



Yamashi 



Keenan 



Yamashital (11972 



993) 



with 4 and 



19751 ). 1 of 



1972 



1975) th at ha d 16 stars in common with our 



Barnbaum et al. 



Keenanl ( 



1996) with 11. 



Barnbaum et al. 



n total, 18 stars of 



(119961 ) could be used 



19931 ) and 3 of [ 

in comparison with the j-indices obtained through the NSCC methodology. These three 
sc hemes were also u s ed to comp are the MS-ind ex published w ith ours, however , only 6 stars 



of 



Yamashital (11972 



19751 ). 1 of 



Keenan 



(119931 ) and 3 stars of 



Barnbaum et al 



(119961 ) could 



be used. The consistence of the NSCC indice s with the ones published by these authors 
can be seen in Table [7J The 



Barnbaum et al. 



's classification seems to better agree with the 



NSCC optical indices. 

Concerning the circumstellar envelope indices, only the temperature and the SiC-indices 
could be discussed, because no other scheme of classification had presented an opacity index. 
The temperature index obtained through t he NSCC methodology could be compa red with 



the ones of three schemes of classification: 



Yamashital (11972 



19751 ) with 18 stars, 



Keenan 
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( 119931 ) with 4 and 



Barnbaum et al 



were just compared with the scheme of 



(11996) with 11 ones. O n the other hand, the SiC- index 



Sloan et al 



(Il998l ). For this last comparison, we 



considered the equivalence of the indices as described in section 14.21 and the authors have 
12 stars published in common with the ones in our sample. 



6. Conclusions 

This New Scheme of Classification of C-Rich AGB Stars can not be applied to any 
carbon star. It is not a good application for AGB carbon stars that have a very thick 
envelope, e.g. extreme AGB stars can not be treated this way as their optical spectra are 
highly obscured. Moreover, the methodology employed is not tied to the sample presented, 
it has the flexibility to serve to all these kind of carbon stars by using the given coefficients 
and parameters. The indices correspond to either a direct measurement of the intensities 
and equivalent widths of features observed in low resolution spectra of carbon-rich stars or 
obtained through a radiative transfer model fit to infrared data. 

The seven indices presented describe in detail the complex scenario of the carbon rich 
stars. It is possible just by analyzing the compact final notation to get a full set of basic 
information about an AGB carbon star. As all calibrations were established based on the 
well quoted works, the indices and levels employed represent the more successful historical 
parameters on the study of carbon stars. 

Regarding the use of the indices, we demonstrated how the SiC/A.C. ratio and the 
r parameters could together provide information about the evolutionary stages of the 
carbon-rich stars. The evolutionary sequences here proposed may provide a new insight 
over the nature of objects such as the odd C-J stars. 
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Fig. 1. — On the left is the pseudo-continuum used as reference for C2- and j-index param- 
eters, set between maxima A 5722A, A 6202A and A 6620A. On the right is the one defined 
between maxima A 4962A and A 5030A used in the measure of the MS-index parameter. 



- 36 - 




6160 6170 6180 6190 6200 5700 5800 5900 6000 6100 6200 

Wavelength (A) Wavelength (A) 



Fig. 2. — On the left is normal and isotopic C 2 band, A 6192A and A 6168A respectively 
from MC79 2-11, and on the right is a set of C2 absorptions at A 5722A - A 6202A from 
C* 1130. 
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Fig. 3. — The isotopic and normal equivalent width of the CN bands for C* 749, which 
head-bands are located at 13 C 14 N A 6260A and 12 C 14 N A 6206A. 
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Fig. 4.— The MS-band structures in the A 4600A - A 5100A spectral range for W Pic. 
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Fig. 5. — On the left, the Ball and Sri lines are stronger as the CN band in RV Cen spectra, 
but the same does not happens for C* 1130. On the right, the P-branch at A 4352A can be 
seen just in C* 849. 
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Fig. 6. — This figure presents best fit to the models. The open circles denote photometric 
data taken from SIMBAD, full line are ISO data (a) and LRS IRAS data (b), dashed lines 
represent the best models, (a) is the best fit to V636 Mon which presents a 11.3 fxm feature, 
due to silicon carbide (SiC) grains; (b) is the best fit to BM Gem, which presents silicate 
emission features at 9.8 /iin and 18 /xm. At the top right of each graphic it is possible to see 
the emission features in detail. 
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Fig. 7. — On the left, the spread of the sample along the domains of the C 2 -index, while 
on the right it is along the j-index domains. The standards used to determinate the ranges 
of both index are marked with open circles, and the boundary edges are represented with 
a dashed lines. The full line represents the situation in which the primary and secondary 
parameters match. On both figures the targets with higher values of the indices are named. 
On the right it is possible to see a correlation between the parameters adopted for the J-stars, 
although the same is not seen for the ones with j ^ 4. 
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Fig. 8. — C-N stars in a decreasing sequence of their temperature indices. 
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Fig. El — continued 
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El — continued 
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Fig. 9. — C-N stars in a decreasing sequence of their temperature indices. 
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Fig. [91 — continued 
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Fig. [91 — continued 
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Fig. 10. — C-N stars in a decreasing sequence of their temperature indices. 




Fig. [TOl — continued 
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[TUl — continued 
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Fig. 11.- 



C-N stars in a decreasing sequence of their temperature indices. 
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Fig. [TTJ — continued 
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Fig. 12. — C-N stars in a decreasing sequence of their temperature indices. 
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Fig. [12j — continued 
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Fig. [12j — continued 
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Fig. 13.- 



xny pazT{Biujoj\[ 
C-R stars in a decreasing sequence of their temperature indices. 
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Fig. [131 — continued 
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fT3l — continued 
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Fig. 14. — C-J stars in a decreasing sequence of their temperature indices. 
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Fig. [TU — continued 
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Fig. [TU — continued 
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Fig. 15.- 



C-J stars in a decreasing sequence of their temperature indices. 
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Fig. [15j — continued 
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[T5l — continued 
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Fig. 16 



— Best fit model to the 11.3 jim. feature for the sample 
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Fig. HEl- 



continued 
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Fig. [161 — continued 
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Fig. HEl- 



continued 
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Fig. 17 



. — Best fit model to the 9.8 /zm feature for the sample 
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Fig. 18. — The SiC/A.C. x r diagram shows the both parameters calculated through the 
radiative transfer model for the C-N (close circles), C-R (open triangles) and C-J stars (open 
squares) of the sample. Four evolutionary sequences are also represented by exponential 
functions with (3 = 0.02, 0.07, 0.13 and 0.24 associated with SiC / A.C . initial = 0.04, 0.14, 0.26 
and 0.48, respectively. 



Table 1. Ranges of the average indicator for each level of the C 2 -index. 



C 2 X p min <p<p 

max 

C 2 X Pmin < P < Pmax 

X = l 0<p<0.15 ••• X = 5 0.75 <p< 0.95 

X = 2 0.15 <p< 0.35 ••• X = 6 0.95 < p < 1.15 

X = 3 0.35 <p< 0.55 ••• X = 7 1.15 < p < 1.35 

X = 4 0.55 <p< 0.75 ••• X = 8 1.35 < p < 1.55 



Table 2. Ranges of the average indicator for each level of the j-index. 



jx 




Qmin ^ Q ^ Qmax 


jx 




Qmin ^ Q *^ Qmax 


X = 





0<g<0.15 • 


• X = 


4 


0.45 < g < 0.55 


X = 


1 


0.15 <g< 0.25 • 


■ X = 


5 


0.55 < g < 0.65 


X = 


2 


0.25 < g < 0.35 • 


■ X = 


6 


0.65 < g < 0.75 


X = 


3 


0.35 < q < 0.45 • 


■ X = 


7 


0.75 < q < 0.85 



Table 3. Combined inspection for the spectral subclasses. 





Strong 


Weak/ Absent 


Unclear 


Ball A 4554,4934A; Sri A 4607A 


C-N or C-H 


C-R 


C-R or C-N 


P-branch A 4352A 


C-R or C-H 


C-N 


C-R or C-N 


Blue- Violet Extention 


C-R or C-H 


C-N 





Table 4. Ranges of each l evel of the r-index. The la s t colu mn presents the groups from 

Lorenz-Martins fc Lefevrd (I1994I ). 



tX r min < t < T max Opacity Level 



r = r < 0.10 Optically Thin Envelope (I) 

r=l 0.10<r<0.20 

r = 2 0.20 < r < 3.00 Intermediate Opacity (II) 

T = 3 r > 3.00 Optically Thick Envelope (III) 



Table 5. Ranges of each level of the SiC-index. 



SiC X SiC m i n < SiC < SiC max SiC X SiC m i n < SiC < SiC max 



X=l SiC < 0.10 ■•■ X = 3 0.20 < SiC < 0.30 

X = 2 0.10 <SiC< 0.20 ••• X = 4 SiC > 0.30 



Table 6. Ranges of each level of the Temperature- index based on the effective 

temperature of the central star. 



T X Train ^ T <C T max T X Train <^ T <C T, 



max 



X = 3 2Q00K < T < 2800fsT • • • X = 5 2200X < T < 2A00K 
X = 4 2400fsT < T < 2600fsT • ■ • X = 6 T < 2200X 



Table 7. Comparison with Other Schemes 



fYamashita 1972 . 19751 CKeenan 19931 fBarnbaum et al. 1996 ) fSloan et al. 19981 



Spectral subclasses 




2(4) 


11 (11) 




C2-index 


10 (17) 


2(4) 


11 (11) 




j -index 


12 (18) 


1 (1) 


2(3) 




MS-index 


4(6) 


1(1) 


2(3) 




r-index 










SiC- index 








8(12) 


Temperature index 


12 (18) 


1 (2) 


6 (10) 





number of stars with matching classification, followed by the total amount of stars in common with our sample, in 
parenthesis, for each published scheme 



Table 8. Results: New Scheme of Classification of C-rich Stars 



Target Target 



W Pic 


C-J5 


C 2 6 


j7 


MS3 


rO 


SiCl • 


• • C* 1467 


C-J:5 


C 2 4 


j4 


MSI 


rl 


SiC3 


Y Tau 


C-N3 


C 2 4 




MSI 


rl 


SiC2 • 


• • XZ Vel 


C-N:5 


C 2 5 




MS3 


rl 


SiC3 


AA Lep 


C-N:4 


C 2 4 






r2 


SiC2 ■ 


• • C* 1705 


C-H:4 


C 2 4 




MS3 


rO 


SiC4 


BN Mon 


C- 4 








rl 


SiC3 • 


• • U Ant 


C-N3 


C 2 5 




MSI 


rO 


SiC2 


BL Ori 


C-N3 


C 2 3 






rO 


SiC2 • 


• • U Hya 


C-N3 


C 2 4 




MS3 


r2 


SiCl 


V636 Mon 


C-N5 


C 2 5 




MS3 


r2 


SiCl • 


• • TV Car 


C-N 






MS3 


r2 


SiC2 


DH Gem 


C-R2 


C 2 4 






r2: 


SiC4: • 


• • TZ Car 


C-J2 


C 2 6 


j6 


MS3 


r2: 


SiC2 


CR Gem 


C-N5 


C 2 3 




MSI 


rl 


SiC2 ■ 


• • V Hya | 


C-N:6 


C 2 5 




MS3 


r2 


SiCl 


GP CMa 


C-N6 


C 2 2 




MS3 


r2 


SiCl • 


• • S Cen 


C-J2 


C 2 6 


j7 


MS3 


rl: 


SiC4 


NP Pup 


C-R3 


C 2 3 






rO 


SiCl • 


• • SS Vir 


C-N:3 


C 2 5 




MS3 


r2 


SiCl 


CL Mon 


C-J:6 


C 2 4 


j4 




t2 


SiCl • 


• • C* 2031 


C-N:6 


C 2 5 




MS3 


r2 


SiC2 


HV CMa 


C- 6 








r3 


SiCl • 


• • RX Cru 


C-N:6 


C 2 5 




MS3 


rl 


SiC4 


RY Mon f 


C-N5 


C 2 5 




MS3 


rO 


SiC2 • 


• • HD 113801 


C-R: 


C 2 2 




MSI 






W CMa | 


C-N:3 


C 2 4 






rO 


SiC3 ■ 


• • T Mus 


C-J3 


C 2 7 


j7 


MS3 


r2 


SiC2 


BO CMa 


C-N5 


C 2 5 




MS2 


rl 


SiC4 • 


• • V971 Cen 


C-J: 


C 2 5 


j4 


MS3 







Table 8 — Continued 



Target Target 



BM Gem 


C-J4 


C 2 5 


j5 




r2 


Jpec 


• • RV Cen 


C-N3 


C 2 4 




MSI 


rl 


SiCl 


BE CMa | 


C-J2 


C 2 5 


j5 


MSI 


rO: 


SiC4: • 


• • C* 2208 


C-J: 


C 2 6 


j7 


MS3 






C* 749 


C-J5 


C 2 5 


j6 


MS2 


r2: 


Jpec 


• • V Oph 


C-N:4 


C 2 4 




MSI 


rO 


SiC4 


C* 849 


C-H:6 


C 2 4 






r2 


SiC2 • 


• • SU Sco 


C-J5 


C 2 5 


j4 


MS3 


rO 


SiC2 


C* 846 


C-N5 


C 2 5 




MS3 


rl 


SiC4 • 


• • V2309 Oph 


C-N:5 


C 2 5 




MS3 


rl 


SiC2 


C* 1003 


C-J4 


C 2 5 


j4 




r3 


Jpec 


• • TW Oph f 


C-N:5 


C 2 5 




MS3 


rO 


SiCl 


FK Pup 


C-N5 


C 2 4 




MS3 


rl 


SiC2 • 


• • FO Ser 


C-J3 


C 2 6 


j7 


MS3 


r2: 


SiCl 


RY Hya f 


C-N4 


C 2 4 




MS3 


r2 


SiC2 • 


■ • C* 2535 


C-N:6 


C 2 4 




MS3 


r2 


SiCl 


C* 1130 


C-J:4 


C 2 5 


j4 




r2: 


Jpec 


• • UW Sgr 


C-R:* 


C 2 3 










T Cnc 


C-N5 


C 2 5 




MS3 


rl 


SiCl • 


• • BI Cap 


C-R:* 


C 2 4 










MC79 2-11 


C-J6 


C 2 6 


j6 




r2: 


JPec 

















f representative standards of the NSCC 
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